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How a spin-glass remembers. Memory and rejuvenation from intermittency data: an 

analysis of temperature shifts. 
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The memory and rejuvenation properties of intermittent heat transport are explored theoretically 
and by numerical simulation for Ising spin glasses with short-ranged interactions. The theoretical 
part develops a picture of non-equilibrium glassy dynamics recently introduced by the authors, 
which links irreversible 'intermittent' events, or 'quakes' to thermal fluctuations of record magnitude. 
The pivotal idea is that the largest energy barrier &(f„,r) surmounted prior to t^ at temperature 
T determines the rate r^ oc l/t^ of the intermittent events near t^. The idea implies that the 
same rate after a negative temperature shift should be given by r, oc l/t„ '. The 'effective age' 
tw ^ iw has an algebraic dependence on t^ , whose exponent contains the temperatures before and 
after the shift. This analytical expression is confirmed by numerical simulations. Marginal stability 
suggests an asymmetry between cooling and heating, i.e. a positive temperature shift T ^r T' should 
erase the memory b{tm,T). This is confirmed by the simulations, which show a rate Vq oc l/im, 
controlled by the barrier b{tw,T') > b{tw,T). Additional 'rejuvenation' effects are also identified in 
the intermittency data for shifts of both signs. 

PACS numbers: 65.60.-|-a;05.40.-a;75.10.Nr 



I. INTRODUCTION 



Many fascinating aspects of non-equilibrium glassy 
dynamics are revealed by temperature changes applied 
while the system ages, both experimentally, see e.g. |l|, 
Hi Hi 111 and in numerical simulations [1, 0, H- Spin- 
glass experiments often measure the out of phase AC 
susceptibility |3|x"(a;,iu,), or the (related) derivative of 
the ZFC magnetization with respect to the logarithm of 
time S{t,tu]). In the protocol introduced in ref. [3J the 
temperature is decreased at a constant rate, except for 
one temperature To which is maintained for e.g. hours. 
During this aging stage, x" decreases, revealing the char- 
acteristic dynamical stiffening generally present in glassy 
systems. As cooling is resumed, the susceptibility returns 
to the T dependence of an un-aged system. The time 
spent aging at one temperature is thus of little avail at 
a lower temperature, and the x" versus T curve features 
a dip centered at Tq. That this dip is reproduced when 
re-heating at a constant rate through Tq demonstrates 
the persistence of the memory of the configurations vis- 
ited. A later extension of this work ^j shows that two 
dips can be produced and retrieved by aging at differ- 
ent temperatures. More importantly, the memory can be 
erased by overheating, i.e. in a second coohng phases the 
dip is no longer present, provided that the intervening 
temperature has been sufficiently high. 

In temperature shifts experiments, the system is aged 
for a time t^, and a small magnetic field is then turned 
on to probe the dynamics. Under isothermal conditions, 
S{t.tw) has a maximum ai t — tw If a small tempera- 
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ture shift of either sign is imposed at tw , this maximum 
moves jlj. The position of the maximum on the t axis 
defines an effective age, which is studied as a function 
of various parameters. Ref. |j] is particularly concerned 
with the issue of the invertibility of the relation between 
age and effective age, an issue which we return to in the 
discussion. Numerical simulations are also able to follow 
the time dependence of the energy density e{t) after a T 
shift: the so called 'Kovacs' effect is best observed when 
shifting the temperature from Ti to T2 at a late time, for 
which the energy is already close to its equilibrium value 
at T2: After the shift, e{t) approaches the isothermal en- 
ergy density at T2 in a non-monotonic way. This effect 
can be used to define an effective age |3, Q as a time shift 
allowing one to superimpose the the T shifted e{t) onto 
its isothermal counterpart at the final temperature. 

The present investigations give only indirect insight on 
configurational changes and how these changes are mem- 
orized by the aging spin-glass. They do, however directly 
shed light on a related issue, i.e. how the size of the dy- 
namical barriers surmounted in an aging system is mem- 
orized. The data allow a parameter-free theoretical in- 
terpretation by a recent model of non-equilibrium relax- 
ation 9] , which provides a unified description of pure ag- 
ing by attributing a pivotal role to energy fluctuations of 
record size. The simulations combine temperature shift 
techniques with a statistical analysis of heat transfer ^ 
motivated by the intermittency studies of refs. 11, 13.ll3j| 
in soft-condensed matter systems and glasses. In this 
way, many important aspects of aging dynamics can be 
extracted from a one-time quantity, the energy. This is 
potentially useful whenever eliciting a linear response is 
difficult or impossible. The definitions used for the effec- 
tive age are based on the aging behavior of the rate of 
intermittent events, and differ from the definitions used 
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FIG. 1: Left figure (squares): tlie PDF's of the energy 
differences H{St,t,t^) = E{t^ + (n -f l)5t) - E{t^ + nSt), 
n ^0,1, 2... , nSt <t for isotfiermal aging at T = 0.20. Tfie 
three data sets are taken in the observation windows t™ ,tw+t, 
with i™ = 4400, 8800 and 17600. The observation time and 
sampling interval are t = 1200 and St = 22. The PDF is seen 
to approach a Gaussian as the age increases while the other 
parameters remain fixed. Right figure(circles): all parame- 
ters are as above, but the system is now initially quenched 
to T = 0.25 and further cooled to T = 0.20 at f„ = 2140. 
Insert: the skewness of the PDFs as a function of the age at 
which the data collection commences. Squares and circles cor- 
respond to isothermal and temperature shift conditions. In 
the latter case the skewness has a considerably smaller varia- 
tion with the age, i.e. the PDFs differ less among themselves 
than the isothermal ones. The numerical values are smaller 
and similar to those of older isothermal system. 



in studies of linear response effects and energy density, as 
discussed in detail at the end of the paper. Nevertheless, 
whenever a comparison is possible, concurring physical 
conclusion can be drawn. 

The Probability Density Function (PDF) of the heat 
flow data shows that two types of events occur: reversible 
fluctuations, which are described by the Gaussian part of 
the PDF, and large, irreversible events, dubbed quakes, 
through which the excess energy initially trapped in the 
sample is transferred to the heat bath. These events ap- 
pear as an asymmetric exponential tail of the PDF. The 
general form of the PDF is easily understood: thermal 
fluctuations occur independently and simultaneously in 
different part of the systems, whence their sum is Gaus- 
sianly distributed. Occasionally, a thermal fluctuation 
induces a large configurational rearrangement associated 
with a correspondingly large energy release. Such a quake 
is arguably composed of a series of consecutive 'downhill' 
flips of contiguous spins, and its size is exponentially dis- 
tributed if each movement can be the last one with equal 
probability. 

A closer inspection of the data reveals that the rate 
rJtw) of the quakes: i) is temperature independent [ifl 
Il4| . which is striking for activated dynamics, and ii): de- 
creases, under isothermal conditions, as the reciprocal of 
the age 1/tw The effective age fj/^ of a sample, possibly 
subject to T shifts, can hence be defined as l/r^, and its 



behavior can be empirically investigated under different 
conditions. The theoretical idea used to explain these 
findings is that a glassy system 'remembers' the largest 
or 'record' energy barrier b{tw , T) surmounted during its 
(iso)thermal history up to time tw |3;ll3- Specifically, 
only fluctuations exceeding b are able to produce the in- 
termittent events on each time scale. For isothermal ag- 
ing, record dynamics explains the decelerating rate of 
events rg(tu,) ex 1/t^ and accounts for the temperature 
independence of r^. Within this framework, a mathemat- 
ically consistent explanation is only possible if the dy- 
namically available configurations are distributed in en- 
ergy according to an exponential 'local density of states' 
(LDOS), a form which is widely observed in complex sys- 
tems. 

In this paper, we draw a further consequence of these 
ideas, by predicting the effect of a temperature shift 
imposed at t^ on the rate rq{tw,T^T') of intermittent 
events. Unlike the isothermal case, this rate could de- 
pend, beside t^ itself, on the temperatures T and T' be- 
fore and after the shift. For negative temperature shifts, 
the form of i^^-* is derived analytically under the assump- 
tion that the memory of the extremal barrier 6(tu,, T) es- 
tablished prior to the shift persists and that its value still 
determines the evolution at the lower temperature T' . 

A positive shift allows stronger energy fluctuations 
which very quickly overcome 6(t^,, T) and reach b{tw, T'), 
the extremal barrier for isothermal relaxation at the final 
temperature T'. Hence, the intermittent energy release 
continues at a rate r^ ex l/i^,. These properties are fully 
confirmed by the simulations. 

A second observation is that the extremal barrier does 
not account for the proportionality constant between r^ 
and 1/t^-^ (or l/t^)- This constant is larger (i.e. stronger 
intermittency) after shifts of either sign than for the 
isothermal relaxation. In this sense, the shifts make the 
system appear younger, and induce what could be called 
a rejuvenation effect. 

In the next section we summarize the theoretical ideas 
linking non-equilibrium drift to record fluctuations. In 
section IIIII we derive the formula for the effective age, 
which is tested in Section llVl where the numerical results 
are also given. Section^provides technical details on the 
simulations and on the data handling, and section lvTI is a 
brief summary followed by a discussion of how the present 
results relate to previous investigations. 



II. WHY RECORD FLUCTUATIONS COULD 
BE IMPORTANT IN GLASSY SYSTEMS 

Extremal fluctuations are largely irrelevant in equilib- 
rium, but may have pivotal effects in metastable situa- 
tions, i.e. when exceeding a threshold can push the sys- 
tem into a completely different dynamical path. Consider 
flrst, as an example, an artifact, e.g. a skyscraper, which 
collapses if hit by seismic events above a certain magni- 
tude. We assume that, against good engineering practice. 



a 'marginally better' skyscraper is built following any col- 
lapse, which is able to resist events up to the size causing 
the latest break-down, but unable to resist any — even 
slightly — stronger event. In the sequence of gradually 
stronger skyscrapers built according to this prescription, 
each element collapses if hit by an earthquake stronger 
than all its predecessors. The magnitude is irrelevant, 
and the statistics of record-sized earthquakes describes 
in this case physical events. 

Consider now a sequence of random numbers drawn in- 
dependently from the same continuous distribution. The 
distribution of the number n of records, i.e. numbers 
larger than all previous extractions, is independent of the 
distribution from which the data are extracted |l5|. An 
approximate statistical description of the records j^, UM 
is the log-Poisson distribution 



Pn{tl,t2) 



(alog(t2Ai))" 



(iaAi)"", a < 1, t2>h 

. (1) 

where P„ is the probability that n records will occur be- 
tween trials ti, and where the parameter a describes that 
the records can occur in a independent processes run- 
ning 'in parallel'. We shall later think of these trials 
as a time sequence, and of the random numbers as en- 
ergy values produced by thermal fluctuations. Further- 
more we will mainly need the formula in the limit where 
t2 ~ ii. Glossing over the discreteness of these parame- 
ters, and with ^2 = t^ the rate of events near t^, falls off 

as Tq oc l/tw 

Marginal stability provides a physical mechanism by 
which the rank, and not the magnitude of the fluctu- 
ations determines the non-equilibrium dynamics associ- 
ated with changes of attractor. Because the temporal 
statistics of record events is system independent, the con- 
cept is potentially useful whenever metastability is im- 
portant. Its physical relevance was first noticed in a sim- 
ple model of Charge Density Waves dynamics [Ia| , later 
extended to describe the glassy dynamics of various sys - 
tems with multiple metastable configuratio ns JlsL [l7l [l8| . 
and more recently to thermal aging [9l[lOl[l4| of glasses. 
The numerical investigations of ref. 1^ show explicitly 
how selected energy records characterize the attractors 
visited during unperturbed thermal aging after a quench. 

A priori, marginal stability can play a dynamical role 
whenever: i) Numerous attractors, or 'intrinsic states' ex- 
ists , with different degrees of stability, the latter gauged 
by a typical exit time or, equivalently, exit energy bar- 
rier. A marginal increase of the degree of stability of the 
attractor selected can arise from an entropic effect, i.e. 
simply reflect an overwhelming predominance of shallow 
attractors. 

ii) the initial deep quench typically produces configu- 
rations with a large energy excess. Similar to the gravi- 
tational energy of a building released irreversibly during 
its collapse, this energy can be released in relative large 
and (hence) irreversible bursts 9^ JiJj . 

As in the skyscraper example, if the increase in attrac- 
tor stability is marginal, the non-equilibrium dynamics 



is driven by fluctuation extrema. The number of quakes 
occurring in the interval [l,t«,] is then described by the 
log-Poisson distribution given in Eq. ^ with ti ~ 1 and 

To reach quantitative predictions, some basic aspects 
of real-space morphology must also be considered: Aging 
in spin-glasses with short-ranged interactions is charac- 
terized by a slowly growing length scale, the thermal cor- 
relation len gth, which remains shorter than a few lattice 
spacings [a,ll3 for accessible time scales. It is therefore 
reasonable to assume that in an extended system with 
short range interactions many thermalized sub-domains 
exist, whose linear size matches the thermal correlation 
length and which are separated by frozen degrees of free- 
dom. Their reversible thermal fluctuations are statisti- 
cally independent and are described, through the usual 
tools of statistical mechanics, by a 'local' energy land- 
scape comprising, for each domain, the configurations 
available. 

The model assumption explored in the following is that 
an extremal fluctuation within each sub-domain is able 
to trigger a quake 9, 10]. This has two consequences: 
Firstly, the total number of quakes is a sum of indepen- 
dent Poisson processes, and hence itself a Poisson pro- 
cess with average a > 1. This average increases with the 
system size and attains the smallest possible value a — 1 
only if a single sub-domain is present. Secondly, as we ar- 
gue below, combining the thermal nature of the extremal 
fluctuations with the properties of record statistics means 
that the sub-domains must have a nearly exponential lo- 
cal density of states LDOS, of the type observed numer- 
ically in a variety of small glassy systems jl^, I20I l2ll | . 

Ignoring for simplicity any difference between the sub- 
domains, we write this LDOS as 



V{e) ocexp(e/eo). 



(2) 



and note that the condition T < eo is required in order 
to ensure the thermal metastability of any attractor to 
which this LDOS applies. 

Importantly, a log-Poisson description is only possible 
as long as a sharp distinction between local reversible 
fluctuations and irreversible changes is meaningful. The 
limitation does not seem unduly restrictive in deeply 
quenched glassy system, which for many decades of time 
remain far from true thermal equilibrium. The statisti- 
cal independence of fluctuations which is required by the 
theory, is approximately fulfilled at low temperature, as 
trajectories mostly dwell at the bottom of local wells. 

The above ideas can be directly tested in intermittency 
measurements, where different parts of the same data 
set provide the required information |T0|. E.g. the left 
panel of Fig. ^ shows the Probability Density Function 
(PDF) of the heat transferred over a small time inter- 
val St under isothermal conditions. The Gaussian part 
of the PDF describes the reversible fluctuations and is 
found to be independent of tw [lOJ, indicating that the 
attractors successively explored are similar with respect 
to their equilibrium energy fluctuations, as already im- 
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FIG. 2: Main figure, lower data sets: the PDF's of tlie en- 
ergy differences H{5t,t,tw) = E{tw + {n + l)5t) — E{t^ + n5t), 
n = 0, 1,2. . .. Tlie five sets correspond to ages t^ ~ lOOSt 
witli 5t = 10, 20, 40 . . . 160 and observation window i = 1^/2. 
Tfie temperature is T = 0.20. Upper data sets: One ad- 
ditional cooling step takes the system from T = 0.25 to 
T — 0.20. For the i'th data set, i = 1,2, ...5 the change 
occurs at age f™;, determined such that the corresponding 
effective ages are tl^^ = 2000, 4000, 8000, 16000 and 32000. 
The corresponding sampling intervals are St = i„ /lOO and 
t = f5/-'^/2. Keeping the ratio St/f^^ constant collapses the 
data, whence the effective age plays the same role as the ac- 
tual age does in the isothermal case. The insert shows that 
the skewness of the PDF's falls into two distinct groups, with 
only a small variation within each group. 



plied by Eq. |21 even though they may differ with respect 
to their stabihty against fluctuations. The temperature 
dependence of the variance of the Gaussian fluctuations 
agrees vi^ith the equihbrium value calculated from Eq. |21 
lending quantitative support to this assumption 10] . 

The exponential tail describes the intermittent events, 
which are rare but far more frequent than Gaussian fluc- 
tuations of comparable size. The tail has a clear age 
dependence, which is visible in Fig. Q] and which gets 
stronger the smaller t^ is. An increase of the weight 
of the tail is hence tantamount to a rejuvenation effect. 
The rate of intermittent events falls off as r^ oc 1/i^, 
as demonstrated by the collapse of a set of PDF's ob- 
tained for different t^ and 5t with constant ratio btjt^. 
In what follows, we use the collapse of heat transfer PDF 
to draw inferences on the rate of intermittent events, and 
hence on the 'effective age' of the system, after a small 
temperature shift. 



III. HOW A SPIN-GLASS REMEMBERS 

If uncorrelated energy records trigger the quakes, the 
time intervals between these quakes, i.e. the residence 
time tr spent in metastable states, grows in propor- 
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FIG. 3: The figure shows how the quality of the data collapse, 
of the sort shown in Fig. |21 deteriorates as the condition of 
constant St/tll^ is violated. The effect is achieved by chang- 
ing the value of St, which amounts to varying t^' around 
its theoretical value given by Eq. |S] The abscissa is the rel- 
ative deviation, and zero thus corresponds to the theoretical 
prediction. The following negative T jumps are considered 
0.25 -^ 0.20 (circles), 0.35 -^ 0.30 (squares), 0.45 -^ 0.40 (di- 
amonds)and 0.55 — > 0.50 (stars). Lines are guides to the eye. 
Note that in all cases the best collapse is obtained at or very 
near the predicted value of the effective age. 

tion to the age Q and is independent of the tempera- 
ture. This so-called 'pure scaling' behavior is well known 
and analytically predicted for log-Poisson statistics. We 
presently choose to ignore the usually small but always 
systematic deviations from pure aging, sub- or super- 
aging, which appear in both numerical and experimental 
results. 

A second expression for t^ can be derived using that 
the system is thermally equilibrated in the metastable 
states between the quakes. In this case, hopping over the 
extremal barrier b{tw , T) (see below) is mediated by equi- 
librium thermal ffuctuations. Hence, the residence time 
in a metastable state entered on a time scale t^ is propor- 
tional to the Arrhenius scale tr = ta c<. eKpibitu,, T)/T). 
Consistency with our previous argument requires ta ~ tw 
which is mathematically possible only if the LDOS from 
which the energy fluctuations are drawn has the expo- 
nential form given by Eq. |21 

Equation [21 describes the distribution in energy of the 
configurations accessible to an arbitrary thermalized sub- 
domain, with zero defined as the lowest available energy. 
The probability of visiting states of energy larger than e 
without exiting the attractor is then 



Pe9(e) = exp(-ea(r)), 
where the reduced inverse temperature 

a(T) = (l/r-l/eo) 



(3) 



(4) 



describes, for temperatures well below eq, the entropic ef- 
fect of the exponential divergence of the LDOS. Through- 



out, the value cq = 0.86 is used. This value, which is close 
to the critical temperature of the 3d Edwards- Anderson 
spin-glass 22|, was obtained in ref. [i3| by fitting the 
theoretical form predicted by Eq. |31to the heat capacity 
data inferred from the Gaussian heat fluctuations. 

The largest fluctuation occurring during (0, tw] is the 
largest of 0{t^) energy values drawn independently from 
the exponential equilibrium distribution O) . According 
to a standard result of extremal statistics [23 (NB: this 
result only requires that the tail of the Boltzmann dis- 
tribution be exponentially distributed) its average value 
grows with the age as 



b{t^,T) = ln(t^,)/a(r) + const. 



(5) 



whence the Arrhenius time ta — exp(a(r)6(itu, T)) re- 
mains equal to i^,, as anticipated. 

The above argument is now extended to the situation 
where the temperature is instantaneously lowered from T 
to T' at age t^: the extremal barrier 6(tu,, T) previously 
established remains the same, but must now be overcome 
at a lower temperature in order to induce a quake. The 
relevant Arrhenius time is then given by 



TA 



-+eff 



t^Jf cx exp(a(rXi„,T)) oc t 



AT') 



T' < T, (6) 



where Eq.[Slis used for the right-most equality. The rate 
of intermittent events after the shift is Vq oc l/f^^ . Note 
that the argument ignores the decay of the rate occurring 
while the statistics is collected. In the cases considered, 
the drift remains small since the observation time is well 
below the effective age. 

Equation 1^1 implies that the PDF's of the heat trans- 
fer obtained for various 5t and i^, values collapse if the 
ratio St/f^f is kept constant. This is demonstrated in 
Figs.[21and|31 which show that t^f^ plays a role analogous 
to t^, in isothermal aging experiments (See Figs.[21and0] 
and ref. |10J). As anticipated, the step has an additional 
effect: the collapse leads to a master curve with a con- 
siderably stronger exponential tail than the curve corre- 
spondingly obtained from isothermal data. 

In summary, a temperature drop has two effects: a: 
It effectively increases the size of all barriers, leading to 
an effective age greater than the age. In this respect 
the system looks older, b: The amount of heat delivered 
by the intermittent events is higher, whence the system 
looks in this respect younger. As later discussed, a simi- 
lar duality is seen in refs. [5|,|^{. A positive temperature 
shift qualitatively differs from a negative shift. The ex- 
tremal barrier 6(t„,r) still sets the time scale for the 
first intermittent event after the step. Since this barrier 
is now strongly reduced relative to the typical size of the 
thermal fluctuations, its crossing is neither rare nor nec- 
essarily irreversible. The additional time needed to reach 
the extremal barrier b(t.w,T') is, clearly, bounded by t^, 
whence we can expect the data to be scaled by keeping 
St/ttj^: constant. This is confirmed by Fig. 0] showing an 
excellent data collapse after the positive shift. Again, 



comparing the master curve with its isothermal counter- 
part shows a stronger exponential tail. Hence, a positive 
step rejuvenates the dynamics in the same sense as dis- 
cussed in item b above. 



IV. SIMULATION RESULTS 

The simulation results are now described in more de- 
tail. As in ref. |1^, the statistics of the amount of heat 
exchanged over a time interval of length St < t^ is col- 
lected for n contiguous intervals, with n chosen such that 
nSt =^ t. The fc'th value sampled is thus 

Hk^E{t^ + {k+l)St)- E{t^, + kSt), k<n, (7) 

where E is the energy of the system. Since the statis- 
tical properties of the heat transfer depend on St, t and 
t^, only, the process itself is denoted by H{St, t, t^)- The 
drift toward lower energies is described by an asymmet- 
ric exponential tail, and the fiuctuations by the Gaussian 
part of the PDF. For fixed t and St, the PDF's approach 
a Gaussian shape as t^ increases and the skewness (third 
central moment) correspondingly approaches zero from 
below. This is seen in the left frame of Fig. Q] and in 
the corresponding insert. The three PDF's are taken 
at ty, = 4400, 8800 and 17600, all with with t = 1200, 
St = 22 and T = 0.20. The data shown in the right hand 
panel of Fig. ^ are obtained with the same parameters, 
except that the temperature is instantaneously decreased 
from T = 0.25 to T = 0.20 at U = 2140. The Gaussian 
fluctuations are not affected by the step. The tails, how- 
ever, are less pronounced than in the isothermal case, 
and thus qualitatively similar to those of an isothermally 
aged but older system. 

Figure 12 contains twelve data sets, which fall into 
two distinct groups of equal size. Data with the most 
pronounced intermittent tails are taken after instanta- 
neously decreasing the temperature from T = 0.25 to 
0.20 at ages U = 275, 460, 770, 1283, 2143 and 3578. Ac- 
cording to Eq. |S1 these values correspond to the effective 
ages fjf = 2000, 4000,... 64000. The data collection 
starts immediately after the step and continues in each 
case up to the age corresponding to {i/2)t'^^ . The en- 
ergy differences in the r.h.s. of Eq. |7| are taken over a 
time interval S — t^f /lOO. For each set of physical pa- 
rameters, 6000 independent runs were performed. The 
collapse demonstrates that the rate of intermittent events 
is proportional to the theoretical prediction l/tf^^. For 
comparison, six data sets were collected isothermally, 
from age t^, to age (3/2)^^,, with St = t^/WO and 
t^, = 2000,4000,8000, ...64000. They collapse among 
themselves but, even though the t^ values pairwise cor- 
respond to the t^^f^ values for the temperature step, the 
two groups of data remain clearly distinct. This is the 
dual effect previously mentioned. 

Very similar results were obtained for other tempera- 
ture steps , i.e. 0.55 ^ 0.50, 0.45 -^ 0.40 0.35 -^ 0.30. 
Figure 13 shows that the quality of the results is highly 
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FIG. 4: Main figure, upper data sets: the heat trans- 
fer PDF's collected after a temperature step imposed at 
tn, = 1000, 2000, 4000, 8000 and 16000, i.e. in the interval t» to 
(3/2)iiu. The energy differences are taken over 5t = 4^/100. 
Keeping St/t^ constant is seen to achieve an excellent col- 
lapse of the data. Lower data sets: the same quantities are 
obtained under isothermal conditions at T = 0.25. The insert 
shows that the skewness values of the data sets fall into two 
clearly distinct groups. 



sensitive to the value chosen for f^J"^ , with the best results 
obtained at or near the theoretical value. The ordinate 
empirically quantifies the quality of the data collapse of 
two PDF's, say a and b, and is obtained as follows. First, 
the extreme tails, which are mostly affected by statistical 
error, are cut off. We then let p{a) and p{b) denote the 
probabilities carried by the rest of the tails. The error is 
defined as e{a,b) =| {p{a) - p{b)) \ {p{a)p{b))~^^'^. For 
each data set, p is calculated over an interval extend- 
ing over 8 units of energy, i.e. from —4 to —12 for the 
0.25 -^ 0.20 data. This interval is gradually moved to 
the left (more negative values) for 0.35 -^ 0.25 etc., to 
match the increasing width of the Gaussian part of the 
curves as the temperature increases. 

The upper data set in Fig. 0] shows the data collapse 
obtained after a positive temperature shift. As was the 
case for the negative shift, the Gaussian parts are unaf- 
fected, while the intermittent tails are strongly affected. 
The set, comprises six different PDF's from systems of 
different ages, i„ = 2000,4000, . . .32000, with the ratio 
6t/tuj = 1/100. Data are collected from t^ to {3/2)t^, 
and 2000 independent runs were performed for each set of 
physical parameters. Since, as the data collapse demon- 
strates, the rate of intermittent events remains propor- 
tional to 1/iu,, the stronger fluctuations available after 
the shift are able to quickly overcome all barriers in the 
range between &(*„,, T) and 6(itu, T'), the extremal barrier 
which would have been reached isothermally at the higher 
temperature. The latter barrier hence sets the rate for 
the intermittent events. The lower data sets comprise 
five isothermal PDF's obtained under the same condi- 
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FIG. 5: The figure explores how the quality of the PDF 
collapse after a positive temperature step — of the sort shown 
in Fig. |1] — deteriorates as the criterion of constant St/tw is 
systematically violated. The abscissa is the relative deviation 
of the scaling parameter used from tw . The following positive 
T jumps are considered 0.20 — > 0.25 (circles), 0.30 — > 0.35 
(squares), 0.40 — * 0.45 (diamonds) and 0.50 -^ 0.55 (stars). 
Lines are guides to the eye. 



tions, but for fixed temperature T = 0.25. The presence 
of two distinct master curves shows, again, the presence 
of an additional rejuvenation aspect. 

Finally, Fig. O explores how the quality of the collapse 
of PDF's obtained after a positive temperature declines, 
as the criterion St/t^ =constant is systematically vio- 
lated. These results, obtained for initial temperatures 
T = 0.20, 0.30 .. .0.50, upward shifted in each case by 
0.05, show that the conclusion drawn from Fig. ^remains 
valid for a range of temperatures. 



V. TECHNICAL DETAILS 

All simulations are performed for a 16'^ Edwards- 
Anderson spin-glass model with nearest neighbor Gaus- 
sian couplings j24] of unit variance and zero mean. The 
Boltzmann constant is set to one. 

We employ an event driven simulation technique, the 
Waiting Time Method (WTM), also used in the land- 
scape exploration |l4| . The 'intrinsic' or current time 
of the system is initially zero and corresponds for large 
system sizes to the usual Monte Carlo sweep. A slight 
modification of the WTM procedure described in ^^ is 
required to deal with temperature changes beyond the 
initial instantaneous quench at t^ = 0, if the mathemat- 
ical equivalence of the WTM to the standard Metropolis 
algorithm is to be maintained. This is discussed below 
together with the main points of the algorithm. 

Let SEi be the energy difference produced by flipping 
the i'th spin. Initially, random waiting times are drawn 
independently from exponential distributions with aver- 
age Ti = ma.x{exp{5Ei/T),l). These random quantities 



also define, for each spin i, its flipping time. This is 
the time at which the spin would flip, were 5Ei to re- 
main constant in the meanwhile. The simulation iterates 
the following two steps: i) the current time is set to the 
earliest flipping time, and the corresponding spin, e.g. 
number fc, is flipped; ii) the waiting times of spin k and 
its nearest neighbors are re-calculated taking the new 
6Ei/Ts into account. Adding these waiting times to the 
current time yields the updated flipping times of spin k 
and its neighbors. All other spins are unaffected by the 
move and need not be updated, as they are engaged in 
a memory-less (Poisson) waiting process. The modiflca- 
tion alluded to is that, whenever T changes, all 6Ei/T 
values are affected and all the waiting and flipping times 
are re-calculated. 

The data shown are obtained, for each choice of phys- 
ical parameters, from 2000 to 6000 independent runs. 
Each run typically samples 50 values of H (this fig- 
ure varies slightly according to the parameter settings). 
PDFs are thus estimated as the empirical frequencies of 
TV = 1.5 • 10^ data points for 150 bins covering the range 
of the observations, and normalized for convenience to 
cover a unit area. The statistics of the Gaussian part of 
the data is excellent and does not need further discussion. 
The statistical error on the exponential tail is, unfortu- 
nately, not completely negligible. Consider a data bin 
whose theoretical probability is p, with p << 1. The 
number n of independent measurements which out of 
N >> 1 trials fall into this bin is well approximated by 
a Poisson variable, with average and variance both equal 
to pN. The estimate n/N for p has mean p and standard 
deviation (p/iV)^/^. The relative error on a data point 
with theoretical probability p is the ratio of the standard 
deviation to the mean and equals {1/pNy/^. The slow 
convergence with N prevents a perfect collapse of the tail 
data. A relative error of 10%, which is still graphically 
evident in a semi-logarithmic plot, is obtained for our 
range of N for values oi p above 10"'^. 

To skewness (third central moments), which is used 
as a measure of the deviation from the equilibrium-like 
Gaussian shape, amplifies the noise in the poorly sampled 
extreme tail of the PDF. To mitigate the effect, we esti- 
mated the skewness using PDFs curtailed at 10~^. This 
appears to be sufficiently accurate for comparing PDF's 
describing different physical situations. 



VI. SUMMARY AND DISCUSSION 



Statistical analysis of mesoscopic noise is a power- 
ful probe of a-gin g d ynamics in a variety of glassy sys- 
tems [l^ [ll|, U^ [13 1 ss reversible equilibrium fluctua- 
tions and irreversible drift appear as distinct components 
of the signal's Probability Distribution Function (PDF): 
the former are described by a Gaussian of zero average, 
and the latter by an 'intermittent' exponential tail. Our 
investigations utilize the intermittent tail to describe the 
aging behavior of intermittent events, or quakes, showing 



how their rate decays with the age and how the decay is 
modifled by temperature shifts. A unified explanation 
of both these effects and of pure aging behavior |3, |i3 
is achieved by postulating that quakes occur whenever 
a record sized thermal fluctuation crosses the extremal 
barrier established in the course of the preceding ther- 
mal history: a spin glass is thus claimed to remember a 
sequence of ever growing extremal barriers. This infor- 
mation is erased by a positive T shift but preserved by 
a negative shift. The mechanism suggests the presence 
of a strong entropic bias favoring shallow attractors, i.e. 
ensuring that the least stable attractor is selected with 
high probability after each quake. 

The rate of intermittent events following a T shift is 
thus a sensitive probe of the size of the barriers sur- 
mounted by an aging system before the shift, but car- 
ries no information on e.g. the energy released through 
the quakes. The latter crucially affects the energy den- 
sity e(i) as a function of the age, whose behavior under 
T shifts has been thoroughly investigated 0, Q • After 
temperature shifts, e(t) approaches the isothermal en- 
ergy density curve obtained at the final temperature. The 
approach is non-monotonic for positive shifts, a feature 
generally called a Kovacs effect. As the energy density is 
closely related to the heat transport PDF a comparison 
is very instructive. 

The effective ages as presently defined depend on other 
dynamical properties than those extracted from e(t), i.e. 
on the rate of intermittent events rather than the full 
shape of the energy decay. In particular, they do not 
fully describe the dynamical effect of the T shift, e.g. 
the shape of the collapsed PDFs differs for T shifted and 
isothermally aged systems. By contrast, in energy den- 
sity studies the effective age is constructed to make the 
T-shifted and isothermally aged curves coincide. Fur- 
thermore, our definition most directly relates to the bar- 
riers scaled prior to the step and to the effect of the step 
itself, unlike the effective age extracted from e{t) curves 
which most directly gauges the behavior after the ther- 
mal step. 

In an aging system one expects the configurational 
memory to vanish on a time scale equal to the system age. 
For a negative T step we expect the memory of the step to 
vanish on a time scale of the order of f^-^ >> t^,, which 
qualitatively agrees with the results of refs. ja, Q This 
reasoning does not apply to a positive step: the inter- 
mittency rates immediately adjust to those of an isother- 
mally aged system, while the energy density shoots up 
and then appraoches the isothermal curve on a time scale 
which is shorter than t^ but differs from zero. Clearly, 
an abrupt change in the barrier structure does not log- 
ically imply an equally abrupt change in energy or con- 
figuration, and the data tell us that the two types of 
change happen on different scales. In summary, while a 
naive correspondence cannot be established, a substantial 
physical similarity emerges when notational and method- 
ological differences are properly accounted for. 

The same applies to the effect of T shifts other phys- 
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ical properties, i.e. the AC susceptibility x of Edwards- 
Anderson spin-glasses in four and three spatial dimen- 
sions studied in simulations by Berthicr and Bouchaud f| 
and Takayama and Hukushima 6J. These authors pa- 
rameterize their results with effective ages extracted in a 
way similar to the energy density analysis just discussed. 
As in our findings, negative shift have a double effect: 
the apparent rejuvenation consisting in a fast decay of 
x" immediately after the shift and resembling the evo- 
lution right after the initial quench. Secondly, the curve 
'undershoots' the isothermal data at the final temper- 
ature, which makes the system look older. In ref. [g, 
the relation between age and effective age is nearly a 
power-law, i.e. similar to Eq.|Hl The exponent is however 
T/T, and not our a{T')/a{T), where a{T) is given in 
Eq. ^ A positive shift is there described by an effective 
age shorter than the age, while our effective age simply 
remains equal to i^j- This difference corresponds to the 
difference already discussed in connection with the e(t) 
data, but the asymmetry between positive and negative 
shift found in ref. jg is similar to our findings. 

Comparing simulation studies with experimental re- 
sults, as we shall do next, requires an element of faith, 
particularly so for dynamical effects: simulation algo- 
rithms such as the Metropolis algorithm and the WTM 
presently utilized are designed to sample equilibrium dis- 
tributions, but the time scales of e.g. partial equilibra- 
tions are algorithm dependent and cannot be meaning- 
fully connected to experimental time scales. A large em- 
pirical evidence nevertheless points to the relevance for 
dynamical phenomena of the information gained from 
simulations of glassy systems, perhaps unsurprisingly, 
given that the dynamical features experimentally defin- 
ing glassy behavior are similar across a wide range of 
microscopically very different systems, and that aging 
dynamics is, broadly speaking, controlled by time scales 
ratios, rather than magnitudes. 

The temperature and age dependence studies of the 
imaginary part of the AC susceptibility x" |^ i show how 
the memory of configurations explored at a certain tem- 
perature survives cooling and is retrieved by re-heating. 
This information is not directly available from heat trans- 
fer intermittency data. However, the related property 
that the memory is erased by a sufficiently strong over- 
heating [g is qualitatively similar to the asymmetry be- 
tween positive and negative step already discussed. 

The relationship between effective and actual age un- 
der positive and negative temperature shifts is investi- 
gated experimentally in ref. W\. For sufficiently large T 
shifts, a lack of reversibility of the T step, called a 'non- 
accumulative aging effect' is present. This is qualitatively 
similar to our findings, where the shifts applied are of the 
order of 5 % of T^, i.e. large compared to those used in 
ref. |j]. We can thus confirm our conclusion that the 
physical insight obtained from the intermittency anal- 
ysis broadly concurs with the results of other types of 
investigation. 

As intermittent heat transfer says only little about con- 



figuration changes, information regarding the high tem- 
perature sensitivity of the thermal correlation length, so 
called 'chaos' |2a|, which is often linked to rejuvenation 
effects, can only be gained indirectly. We note however 
that since memory and rejuvenation effects appear in sys- 
tems of very modest size, properties of large scale thermal 
objects, e.g. droplet excitations out of a ground state [23 
are largely immaterial for the present discussion. 

As recently shown |lj| , geometrical properties of the 
energy landscape of a spin-glass depend in a systematic 
way on the extremal energy barrier crossed. More gen- 
erally, a linear relationship between the average energy 
barrier crossed and the configurational changes induced 
by the crossing has been measured experimentally [3 
and numerically, in different types of spin-glass mod- 
els [ij,|2^ and using different exploration techniques |29l |. 
Due to the strong connection between barriers overcome 
and configuration changes we can expect that erasing the 
memory of the extremal barrier would also lead to some, 
more gradual, loss of configuration memory. 

A linear relationship between barriers and Hamming 
distances is included in descriptions of spin-glass dynam- 
ics as a master equation on a set of hierarchically orga- 
nized states [^ yH |32, IH, |3J| . The tree models used 
in these descriptions are energy landscape descriptions 
of local regions of a complex system, and are able to re- 
produce salient features of aging, including t/t^ scaling, 
and the effects of temperature shifts based on thermal 
hopping, rejuvenation effects have also been found [33 
in the so-called LS tree l33J, based on the idea that a 
kinetic bias exist favoring shallow attractors. The record 
induced dynamics scenario presently advocated also re- 
quires the presence of a hierarchy of inequivalent dynam- 
ical attractors. However, unlike the case of a static tree 
model, these attractors are not given in advance, but are 
dynamically selected by record energy fluctuations. 

Unlike other methods, the intermittency probe of ag- 
ing dynamics has a direct analytical interpretation, and, 
at least for spin-glasses, appears to offer a different age 
and temperature window where it can conveniently be 
applied. I.e. intermittency rates can be obtained for all 
ages, whether isothermal evolution is fast or slow, but are 
most conspicuous at low temperatures (below T « 0.6 
for the present system) where Gaussian fluctuations are 
smaU [13. 

Secondly, since intermittent events entail large config- 
urational rearrangements, they expectedly influence in 
similar ways different physical properties. E.g. studies 
of the electrical voltage fluctuations in Laponite, a gel, 
and in a polycarbonate glass [33 show that the inten- 
sity of the intermittent signal decreases while the sys- 
tem ages (as in our fig. ^. In spin-glasses, where the 
spin contribution to the energy is not easily separated 
out, one could study the intermittency property of the 
magnetic noise or, in metallic spin-glasses, the electrical 
fluctuations associated with spin flips 36] . The method 
presented can in principle be applied to analyze intermit- 
tency rates of any relevant and experimentally available 



quantity. Measuring the rate of intermittent events as a 
function of the age would be extremely interesting from a 
theoretical point of view, but also extremely challenging. 
Some changes of the PDF's cause by temperature shifts, 
see e.g. Fig. ^ might be observable even in the lack of 
precise rate measurements. 
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